DISTRIBUTION  STATEMENT  A:  Approved  for  public  release;  distribution  is  unlimited. 


Surface  Wave  Processes  on  the  Continental  Shelf  and  Beach 

Thomas  H.  C.  Herbers 
Department  of  Oceanography,  Code  OC/He 
Naval  Postgraduate  School 
Monterey,  California  93943 
Tel:  (831)656-2917 
FAX:  (831)656-2712 
Email:  thherber@nps.edu 

Tim  T.  Janssen 
Department  of  Geosciences 
San  Francisco  State  University 
San  Francisco,  California  94132 
Tel:  (415)338-1209 
FAX:  (415)  338-7705 
Email:  tjanssen@sfsu.edu 

Robert  T.  Guza  and  William  C.  O’Reilly 
Scripps  Institution  of  Oceanography 
Integrative  Oceanography  Division 
La  Jolla,  California  92093 

Tel:  (858)  534-0585  (Guza),  (858)  534-6258  (O’Reilly) 

FAX:  (858)  534-0300 

Email:  rguza@ucsd.edu;  woreilly@ucsd.edu 

Award  Numbers:  N0001409WR20029,  N0001409WR20180,  N00014-07-1-0365, 
N000 1 4-07- 1  -0402,  N000 14-09- 1  -0363 
http://www.oc.nps.navy.mil/wavelab/ 


LONG-TERM  GOALS 

Wind  waves  and  swell  dominate  the  hydrodynamic  and  sediment  transport  processes  on  many 
continental  shelves  and  beaches,  drive  the  nearshore  circulation  and  morphology,  and  play  an 
important  role  in  remote  sensing  applications.  Accurate  wave  prediction  in  coastal  environments  is  a 
daunting  challenge  because  waves  are  affected  by  many  physical  processes,  including  scattering  by 
seafloor  topography,  wave-current  interactions,  nonlinear  effects,  wave  breaking,  and  friction  in  the 
bottom  boundary  layer.  Several  of  these  processes  are  poorly  understood  and  existing  wave  prediction 
models  rely  on  parameterizations  and  empirical  calibration  to  represent  them.  The  long-term  goals  of 
this  research  are  to  obtain  a  better  understanding  of  the  physical  processes  that  affect  ocean  surface 
waves  in  the  coastal  environment  and  develop  improved  wave  prediction  capability. 
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OBJECTIVES 


•  Observe  and  predict  the  seafloor  damping  effects  on  ocean  surface  waves  in  sandy  and  muddy 
coastal  environments. 

•  Advance  deterministic  and  stochastic  modeling  capability  for  nonlinear  wave  evolution  in  coastal 
regions  with  complex  seafloor  topography  and/or  strong  currents. 

•  Improve  the  representation  of  source  terms  in  operational  wave  prediction  models. 

•  Predict  the  nonlinear  shoaling  transformation  of  waves  on  beaches  including  the  excitation  of 
infragravity  motions. 

APPROACH 

By  combining  theoretical  advances  with  numerical  models  and  field  observations,  we  investigate  the 
physical  processes  that  affect  ocean  surface  waves  on  continental  shelves  and  beaches.  The 
transformation  of  wave  spectra  is  predicted  with  models  that  include  the  effects  of  refraction, 
scattering  by  wave- wave  and  wave-bottom  interactions,  and  parameterizations  of  bottom  friction,  and 
wave  breaking.  Extensive  field  data  sets  were  collected  in  ONR  experiments  off  North  Carolina 
(DUCK94,  SandyDuck,  SHOWEX),  California  (NCEX),  and  the  Florida  Gulf  coast  (SAX04/Ripples) 
to  test  these  models  in  a  range  of  coastal  environments.  Recently,  we  conducted  additional 
experiments  on  the  sandy  Martha’s  Vineyard  shelf  and  the  muddy  Louisiana  shelf.  Analysis 
techniques  applied  to  these  measurements  include  various  inverse  methods  to  extract  directional  and 
wavenumber  properties  from  array  cross-spectra,  higher-order  spectral  analysis  to  detect  nonlinear 
coupling,  as  well  as  standard  statistical  methods  to  determine  empirical  relationships  between  observed 
variables.  The  modeling  efforts  include  deterministic  and  stochastic  models  that  incorporate  quadratic 
and  cubic  nonlinearity  and  wide  angle  diffraction  effects,  suitable  for  application  to  energetic  wave 
environments  with  complex  seafloor  topography. 

WORK  COMPLETED 

During  FY09  we  continued  the  analysis  of  field  observations  of  ocean  surface  wave  attenuation  across 
the  wide  and  muddy  continental  shelf  of  Western  Louisiana.  We  deployed  an  extensive  array  of 
instruments  (Figure  1)  during  February-March,  2008,  that  included  two  directional  waverider  buoys, 
six  bottom  tripods  equipped  with  a  pressure-velocity  sensor  and  a  current  profiler,  and  six  bottom 
tripods  equipped  with  a  pressure  sensor.  The  two-dimensional  array  consisted  of  two  cross-shore 
transects  and  an  alongshore  transect  spanning  a  40  by  25  km  area  in  depths  ranging  from  13  to  4  m. 
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Figure  1.  Field  site  and  array  plan  of  the  Louisiana  Mud  Experiment.  Mud  deposits  from  the 
Atchafalaya  River  are  evident  in  the  satellite  image.  White  lines  indicate  depth  contours  in  m. 


The  dataset  includes  numerous  local  wind  sea  events  with  wave  directions  predominantly  from  the 
south  (i.e.  onshore  propagation)  but  also  fetch-limited  wave  growth  conditions  when  the  wind  was 
directed  offshore  from  land.  To  characterize  the  bottom,  and  in  particular  the  surficial  sediment 
properties,  box  cores  were  collected  at  all  NPS  instrument  sites  (Garcia-Garcia  et  al.,  2008)  and 
analyzed  by  the  Johns  Hopkins  University  (JHU)  MURI  team. 
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Figure  2 .  Wave  conditions  along  the  western  transect  during  the  Louisiana  Mud  Experiment :  Blue 
and  red  dots  indicate  estimates  from  the  offshore  directional  waverider  buoy  and  the  inshore  bottom 

pressure  sensor,  respectively. 
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Preliminary  analysis  generally  shows  a  consistent  decay  of  waves  from  the  deeper  to  the  shallower 
instruments  (e.g.  Figures  2  and  3),  similar  to  earlier  observations  (Sheremet  and  Stone,  2003;  Elgar  and 
Raubenheimer,  2008).  The  wave  spectra  evolution  (Figure  3)  shows  strong  decay  (several  orders  of 
magnitude)  of  high-frequency  wind  sea  spectral  levels  and  weaker  decay  at  the  lower  swell 
frequencies.  These  observations  suggest  that  the  dissipation  is  not  the  result  of  a  direct  wave-bottom 
interaction  (which  would  affect  only  longer- wavelength  waves),  but  possibly  the  result  of  heavy 
sediment  suspension  over  the  entire  water  column  affecting  the  hydrodynamics  of  short  wavelength 
waves  (Sheremet  and  Stone,  2003).  The  observed  decay  was  strongest  at  the  shallowest  site  of  the 
central  transect  where  box  cores  showed  fresh  mud  deposits  (Garcia-Garcia  et  al.,  2008)  and  a  visibly 
murky  sea  surface  was  observed  during  the  experiment  cruises. 


Figure  3.  Observations  (left  panel)  of  wave  spectral  evolution  across  the  inner  Louisiana  shelf  on 
18  February  2008  in  fetch-limited  (offshore  wind)  conditions.  The  data  were  collected  on  the  central 
transect  (Figure  1)  with  sensors  located  in  12m  (circle  markers,  black),  10m  (diamond  markers, 
blue)  and,  6m  (square  markers,  red)  water  depth.  Whereas  the  observed  swell  decay  (the  lower- 
frequency  peak)  agrees  fairly  well  with  model  results  (SWAN,  right  panel),  the  observed  higher- 
frequency  wind-sea  evolution  shows  a  dramatic  suppression  of  wave  growth  in  shallow  water  that  is 
not  predicted  by  the  model.  This  comparison  illustrates  the  critical  role  of  wave-mud  interaction 
processes  in  not  only  the  damping  of  swell  but  also  the  wind-sea  development  in  shallow  water. 

RESULTS 

The  focusing  of  waves  on  tidal  currents  and  shoals  can  create  local  'hot  spots'  of  wave  energy,  where 
large,  steep  waves  are  breaking  resulting  in  hazardous  situations,  in  particular  for  small  vessels. 
Although  the  linear  dynamics  of  these  phenomena  are  well  known,  very  little  is  understood  about 
nonlinear  refraction  of  waves  over  strong  shear  currents,  wave  instabilities  in  such  regions,  and  the 
implications  for  wave  statistics.  During  FY09  we  continued  to  develop  a  numerical  model  for  the 
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Figure  4.  Waves  propagating  over  an  opposing  current  jet  (ebb  tidal  flow  scenario)  predicted  by 
nonlinear  wave-current  model  [Janssen  &  Herbers,  2009],  Nonlinear  effects  in  focal  zone  result 
in  frequency-downshift,  non-  Gaussian  statistics,  and  an  enhanced  likelihood  of  extreme  (freak) 

waves  (see  time  series  inset). 

nonlinear  evolution  of  surface  gravity  waves  through  an  inhomogeneous  medium.  The  model 
describes  the  combined  effects  of  wave  refraction  and  diffraction  and  nonlinear  wave-wave 
interactions  through  scattering  source  terms  in  an  angular  spectrum  evolution  equation  [Janssen  et  al., 
2006,  2008;  Janssen  &  Herbers,  2009].  This  general  approach  can  accommodate  realistic  two- 
dimensional  coastal  features  such  as  an  isolated  shoal  or  a  tidal  inlet,  and  the  numerically  efficient 
pseudo-spectral  solution  technique  is  well  suited  to  Monte  Carlo  simulations  of  statistical  properties. 

A  principal  question  we  hope  to  address  in  this  work  is  how  nonlinear  focusing  over  shear  currents  and 
topography  affects  wave  statistics  in  refractive  focal  areas.  Numerical  simulations  of  random  waves 
propagating  over  an  idealized  shoal  or  opposing  current  jet  demonstrate  that  the  abrupt  variation  in 
wave  energy  associated  with  strong  refraction  can  indeed  cause  a  random  wave  field  to  become 
unstable  and  develop  strongly  non-Gaussian  statistics  and  freak  waves  [e.g.  Figure  4;  Janssen  & 
Herbers,  2009], 

IMPACT/APPLICATIONS 

Existing  spectral  wave  prediction  models  contain  crude,  largely  untested  parameterizations  of  seafloor 
damping  effects.  We  have  conducted  new  field  experiments  on  sandy  and  muddy  continental  shelves 
that  provide  the  much  needed  comprehensive  data  sets  for  improving  these  parameterizations.  We  are 
developing  new  deterministic  and  stochastic  wave  evolution  models  that  incorporate  diffraction  and 
quadratic  and  cubic  nonlinear  effects  for  more  accurate  predictions  in  complex  coastal  environments. 
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